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A  numerical  model  of  beach  profile  evolution  is  developed  to  study  the  onshore  and  offshore  sandbar  migration  under 
different  wave  conditions.  The  integrated  model  consists  of  wave  model,  roller  model,  flow  model,  sediment  transport 
model  and  bed  evolution  model.  In  particular,  the  interaction  between  waves  and  wave-induced  undertow  current  is 
considered  in  the  flow  model.  This  is  achieved  by  adding  a  mean  pressure  gradient  term  into  the  firs -order  momentum 
balance  equation  of  flow.  A  simple  empirical  method  is  also  introduced  to  describe  the  variation  in  eddy  viscosity 
during  one  wave  cycle.  A  preliminary  application  of  the  present  model  shows  good  agreements  of  both  onshore  and 
offshore  sandbar  migration  phenomenon  with  the  laboratory  observation.  It  is  found  that  the  offshore  sandbar 
migration  is  dominated  by  undertow  while  wave  affects  the  transport  rate.  For  the  onshore  sandbar  migration, 
although  sediment  transport  is  mainly  driven  by  wave  and  the  undertow  is  weak,  the  transport  pattern  is  considerably 
modified  by  undertow.  It  is  shown  that  the  combined  effects  of  wave  and  undertow  are  important  in  beach  profile 
evolution. 
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INTRODUCTION 

Coastal  beach  profiles  response  actively  to  the  imposed  wave  and  current  conditions. 
Understanding  and  modeling  beach  profile  evolution  under  varying  wave  and  current  conditions  is 
important  in  coastal  engineering.  Under  energetic  wave  climate,  the  foreshore  becomes  flat  and  a 
sandbar  is  often  observed  in  the  offshore  location.  Under  weak  wave  climate,  the  sandbar  may  move 
onshore  and  the  foreshore  becomes  steep. 

Sandbar  is  one  of  the  most  important  geomorphic  features  in  sandy  coasts,  as  it  serves  a  natural 
protection  for  beaches  by  inducing  more  energy  dissipation  of  waves  above  it.  The  sandbar  is  quite 
dynamic.  It  can  move  onshore  or  offshore  and  grow  or  decay,  depending  on  the  site-specific 
hydrodynamic  conditions  and  the  time-scale  of  these  processes.  A  large  amount  of  research  effort  has 
been  dedicated  to  the  formation  and  migration  of  nearshore  sandbars  (e.g.,  Roelvink  and  Stive  1989; 
Gallagher  et  al.  1998;  Hsu  et  al.  2006;  Ruessink  et  al.  2007;  Walstra  et  al.  2012).  While  the  wave- 
induced  undertow  current  moves  sediment  offshore  and  cause  offshore  sandbar  migration  under 
energetic  waves  (Gallagher  et  al.  1998),  the  skewness  of  nearshore  wave  velocities  and  accelerations 
drives  onshore  sediment  transport  and  sandbar  migration  under  moderate  waves  (Hoefel  and  Elgar 
2003).  In  spite  these  mechanisms  have  been  recognized  and  successfully  implemented  in  various 
models  to  simulate  onshore  or  offshore  sandbar  migration,  simulating  both  conditions  together  with  a 
general  model  setting  is  still  a  challenging  problem  (Van  Rijn  et  al.  2003;  Ruessink  et  al.  2007).  This 
may  partially  be  explained  by  some  missing  or  not  well  described  processes  in  the  models.  Recent 
studies  have  investigated  the  relationship  between  breaking  turbulence  and  sediment  concentration  and 
its  implication  to  undertow-driven  sediment  transport  (Aagaard  and  Hughes  2010;  Grasso  et  al.  2012). 
However,  the  interaction  between  the  high-frequency  wave  motion  and  the  low-frequency  undertow 
and  their  combined  effects  on  sediment  transport  are  still  not  well  understood  yet. 

This  paper  presents  the  development  of  a  numerical  model  of  beach  profile  evolution  in  which  the 
wave-undertow  interaction  are  considered.  The  model  was  tested  with  the  onshore  and  offshore 
sandbar  migration  experiments.  The  distinct  features  of  these  two  processes  and  the  effects  of  wave- 
undertow  interaction  are  discussed. 
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MODEL  DESCRIPTION 
Wave  Model 

Wave  characteristics  are  calculated  using  the  wave-energy  balance  equation: 


8x 

where  £w  is  wave  energy,  cg  is  group  velocity,  Z)w  is  the  energy  dissipation  due  to  breaking  and  is 
calculated  according  to  Battjes  and  Janssen  (1978),  as  modified  by  Janssen  and  Battjes  (2007). 


Roller  Model 

The  roller  energy  is  computed  by  (Stive  and  De  Vriend  1994) 
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where  Er  is  roller  energy,  c  is  phase  velocity,  Dr  is  roller  dissipation,  g  is  gravitational  acceleration  and 
jS=0.1  is  the  roller  slope.  The  wave  set-up  fj  is  estimated  with  the  depth-integrated  and  time-averaged 
momentum  equation  including  roller  effects: 
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where  Sxx  is  wave  radiation  stress,  p  is  water  density,  h  is  the  still  water  depth. 

Flow  Model 

The  intra-wave  flow  velocity  is  modeled  based  on  the  first-order  momentum  equation,  written  as 
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where  u  is  the  instantaneous  flow  velocity,  u-,  is  the  wave-induced  free  stream  velocity,  vt  and  v  are 
turbulent  and  kinematic  eddy  viscosity,  respectively,  p  is  the  mean  pressure,  Df  is  the  energy 

dissipation  in  the  bottom  boundary  layer,  rbls  is  the  additional  mean  shear  stress  induced  by  the 
boundary  layer  streaming.  The  first  term  in  the  right  side  of  Eq.  (5)  represents  wave-induced  pressure 
gradient.  The  third  term  represents  the  mean  horizontal  pressure  gradient  related  to  undertow  (Zhang  et 
al.  2011).  The  fourth  term  is  restricted  in  the  bottom  boundary  layer  and  is  expressed  as  (Reniers  et  al. 
2004) 
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where  3  is  the  boundary  layer  thickness. 

The  time-averaged  and  depth-averaged  turbulent  eddy  viscosity  vtm  is  determined  according  to 
Battjes  (1975): 


=  /,«»(-) 
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where/,,  is  a  constant  taken  as  0.025,  Hrms  is  the  root  mean  square  wave  height,  /z,  is  the  elevation  of 
wave  trough  level  from  the  bed. 

In  this  study,  we  assume  that  the  eddy  viscosity  is  linearly  distributed  over  the  depth  and  is  time- 
dependent  during  one  wave  cycle: 
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u*  =  co  s,  cp  u^+sincp 
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where  u,  is  defined  as  the  equivalent  velocity  for  turbulence  generation  and  is  related  to  the 
instantaneous  value  and  acceleration  of  free  streaming  velocity  and  a  phase  shift  angle  (<p).  u*  has  the 
same  form  to  the  bed  shear  stress  parameterization  of  Nielsen  and  Callaghan  (2003)  but  with  different 
physical  meaning.  In  this  study,  <p  is  calibrated  to  be  26°,  which  is  different  from  the  value  proposed  by 
Nielsen  and  Callaghan  (2003).  This  can  be  explained  by  the  fact  that  in  the  present  formulation  <p  not 
only  describes  the  phase  difference  between  free  steaming  and  the  bed-generated  turbulence  but  also 
takes  the  breaking-generated  turbulence  into  account.  It  has  been  found  that  the  model  results  are 
sensitive  to  the  value  of  tp .  Time  series  of  the  free  stream  velocity  is  computed  with  the  functions  of 
Elfrink  et  al.  (2006). 

At  each  cross-shore  location,  the  mean  horizontal  pressure  gradient  ( -1  / p  dP  / dx)  in  Eq.  (5)  is 
not  known  at  first.  This  term  is  not  directly  estimated.  Instead,  it  is  automatically  obtained  using  an 
iterative  algorithm  which  forces  the  undertow  flux  to  satisfy  the  time-averaged  mass  continuity 
equation: 
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where  T  is  wave  period.  At  the  bottom  boundary,  the  no-slip  condition  for  velocity  is  used.  At  the 
wave  trough  level,  the  shear  stress  formula  of  Deigaard  et  al.  (1993)  is  employed. 

Sediment  Transport  Model 

The  time-dependent  sediment  concentration  and  sediment  flux  are  calculated  with  the  sediment 
advection-diffusion  equation: 


dc  dc  d  (  dc  \ 
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(11) 


where  c  is  sediment  concentration,  ws  is  the  settling  velocity  of  suspended  sediment,  es  is  the  sediment 
diffusivity,  which  is  equal  to  vt.  The  pick-up  function  is  used  as  the  bottom  boundary  condition  and  the 
zero-flux  condition  is  used  at  the  top. 

For  bedload  transport,  a  Meyer-Peter  Mueller  formula  (Ribberink  et  al.  1998)  is  used  that 
functions  the  instantaneous  bedload  transport  rate  with  the  Shields  parameter  {0): 

qh=U^(\d\-0.05)L65^(s-l)gdl  (12) 

where  qb  is  the  bedload  transport  rate,  s  is  the  density  ratio  between  sediment  and  water  (s  =  2.65),  dsn 
is  the  median  grain  diameter  of  sediment. 

The  time-averaged  net  sediment  transport  rate  can  be  calculated  by 

g  =  J  +juc  dz^dt  (13) 


where  q  is  the  total  net  sediment  transport  rate. 

Sediment  Transport  Model 

The  bed  evolution  is  computed  by  the  sediment  mass  balance  equation: 


dzb  _  1  dq 

dt  1-  p  dx 


where  zb  is  bed  elevation,  p  =  0.3  is  bed  porosity. 
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RESULTS  AND  DISCUSSION 

A  preliminary  application  is  shown  to  investigate  the  model’s  overall  skills  for  simulating  both  the 
onshore  and  offshore  sandbar  migration  in  a  laboratory  experiment  (Maddux  et  al.  2006;  Guannel  et  al. 
2007).  In  the  experiment,  the  sandbar  was  firstly  moved  offshore  during  3.5  hours  of  energetic  waves 
with  larger  heights  and  shorter  period,  followed  by  an  onshore  migration  of  11.25  hours  under  actions 
of  weak  waves  with  smaller  heights  and  longer  period,  see  Figure  1  for  details. 


Figure  1.  Experimental  set-up. 


Figure  2  shows  the  comparisons  of  wave  height,  mean  water  level  and  undertow  velocity  for  both 
cases.  The  computed  results  agree  well  with  the  measured  data.  It  is  found  that  the  undertow  velocity 
during  offshore  migration  is  underestimated  near  the  bar  crest.  However,  for  the  onshore  migration,  the 
cross-shore  distribution  of  undertow  velocity  is  well  reproduced.  The  undertow  is  considerably  weaker 
during  onshore  migration  than  offshore  migration. 


ifr- 
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Figure  2.  Comparisons  of  wave  height,  set  up  and  undertow  between  prediction  (continuous  lines)  and 
observation  (circles). 


Figure  3  shows  the  comparisons  of  transport  rates  and  beach  profile  evolution  for  both  cases.  The 
calculated  sediment  transport  rate  shows  a  consistent  trend  of  cross-shore  variation  with  the 
measurements.  It  is  found  that  for  both  cases,  bedload  dominates  the  transport  processes  and  the 
maximum  transport  rate  occurs  near  the  bar  crest.  After  calibration,  the  model  reasonably  reproduces 
the  distinct  features  of  on-offshore  sediment  transport  processes  and  the  observed  sandbar  migration 
under  different  wave  and  undertow  conditions,  implying  that  the  present  model  setting  captures  the 
major  physical  mechanisms  in  these  two  cases. 

The  effects  of  wave-undertow  interaction  are  illustrated  in  Figure  4  and  Figure  5.  For  comparison, 
wave-induced  free  stream  velocity  is  set  to  zero  in  the  model  to  consider  undertow  only.  Undertow  is 
neglected  by  eliminating  the  second  and  third  terms  in  the  left  side  of  Eq.  (10).  It  is  shown  in  Figure  4 
that  waves  affect  the  quantity  and  the  local  gradient  of  the  offshore  directed  transport  rate.  In  particular, 
waves  enhance  the  offshore  transport  and  produces  greater  gradients  of  transport  rate  near  the  bar  crest, 
which  is  consistent  with  the  observation.  This  may  arise  from  both  the  enhanced  intra-wave  near-bed 
velocity  and  thus  the  enhanced  bed  shear  stress.  For  the  case  of  onshore  migration.  Figure  5  shows  that 
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although  the  weak  undertow  can  hardly  move  sand  itself,  it  does  considerably  modifies  the  wave- 
induced  transport  pattern.  If  undertow  is  neglected  in  this  case,  the  model  will  significantly  over¬ 
predict  the  onshore  transport  rate  and  fail  to  capture  the  location  of  maximum  transport  rate.  It  is 
summarized  that  the  combined  influence  of  wave  and  undertow  on  sediment  transport  is  not  a  simple 
summation  of  the  individual  effects,  but  is  rather  a  result  of  the  complex  wave-undertow  interaction. 
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Figure  3.  Comparisons  of  transport  rates  and  beach  profile  evolution  between  prediction  (continuous  lines) 
and  observation  (circles  and  dot  lines). 
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Figure  4.  Comparisons  of  net  transport  rate  including  wave-undertow  interaction  (blue  line)  and  undertow 
only  (green  line)  for  the  case  of  offshore  migration.  The  measured  transport  rate  is  plotted  with  red  dot  line. 


Figure  5.  Comparisons  of  net  transport  rate  including  wave-undertow  interaction  (blue  line)  and  wave  only 
(green  line)  for  the  case  of  onshore  migration.  The  measured  transport  rate  is  plotted  with  red  dot  line. 


CONCLUSION 

A  model  of  beach  profile  evolution  is  developed,  in  which  the  wave-undertow  interaction  is 
included.  The  model  well  reproduces  the  observed  trend  of  on-offshore  sandbar  migration  beneath 
different  incident  wave  conditions.  It  is  shown  that  wave-undertow  interaction  plays  a  significant  role 
in  beach  profile  evolution.  The  combined  influence  of  wave  and  undertow  on  sediment  transport  is  not 
a  simple  summation  of  the  individual  effects.  In  the  present  case,  the  offshore  sandbar  migration  is 
dominated  by  undertow  while  wave  affects  the  transport  rate.  For  the  onshore  sandbar  migration, 


6 


COASTAL  ENGINEERING  2012 


although  sediment  transport  is  mainly  driven  by  wave  and  the  weak  undertow  can  hardly  move  sand 
itself,  the  transport  pattern  is  considerably  modified  by  undertow. 
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